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ABSTRACT 



Crystal oscillators are now available with frequency stabilities 
10 11 

of a few parts in 10 or 10 . One application of such stability 

is in navigation systems, where geographic position may be determined 
by phase -comparing a received frequency-stable signal with a local 
frequency-stable signal. Errors may be introduced, however, by the 
effects of static position and of motion on the frequency stability 
of the oscillators employed. A testing program including static posi- 
tion, swinging motion, vibration, and shock tests was devised to 
determine the extent of such effects. The program is described and 
the results are presented in both tabular and graphic form. 
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1 . 



Introduction. 



Due to improvements in crystal-unit and oscillator design and 
construction during the past ten years, it is now not unusual to 
speak of crystal oscillators with the capability of maintaining their 
frequency outputs constant to a few parts in 10^ or 10*\ Oscil- 
lators of this sort have been employed in very low frequency naviga- 
tion systems, where the phenomenon of long-range, phase-stable VLF 
transmissions coupled with frequency-stable signals at the trans- 
mitters and aboard ship, permits a geographic position to be determined 
by a process of phase-comparison. 

The accuracy of such a navigation system depends upon the 
frequency stability of the oscillators. It is recognized that en- 
vironmental conditions such as temperature, vibration, and shock have 
a direct effect on the operation of crystal oscillators and, therefore, 
every effort is generally made to provide the operating conditions 
as specified by the manufacturer of the equipment. 

In addition to the problems of temperature, vibration, and shock 
it was suspected that motion such as that created by the rolling 
and pitching of a ship or buoy might introduce instability. 

It was the purpose of this investigation to determine the effects 
of such motion on the frequency stability of certain crystal oscil- 
lators as well as the effects of static position, vibration, and shock. 



1 



2 . 



Oscillators Tested. 



Tests were conducted with the following oscillators: 

RF Oscillator 0-76A/U 

Western Electric D-175730-L2 Frequency Standard 
Serial No. 43 

Sulzer Laboratories Model D5 Oscillator 
Serial No. 25 

Sulzer Laboratories Model D5 Oscillator 
Serial No. 19 

Sulzer Laboratories Model D5 Oscillator 
Serial No. 22 

Sulzer Laboratories Model 2.5 Frequency Standard 
Serial No. 256 

A brief description of each type of oscillator is given in the 
following paragraphs. Further information may be obtained from the 



The Western Electric frequency standard produces a frequency 
of 100 kilocycles which, according to the manufacturer, can be 

g 

maintained constant to within approximately one part in 10 per day. 

The crystal, a GT-cut quartz plate with a very low temperature co- 
efficient, is maintained at a constant temperature by a two-stage 
oven and a bridge oscillator circuit. The oscillator weighs 90 
pounds; its dimensions are 19" x 20" x 10". 

The Sulzer Laboratories Model D5 oscillator is a solid-state 
unit with a five -megacycle output and a rated stability of five 
parts in 10^ per day. The crystal is an AT-cut quartz plate, vibrating 
in a fifth-overtone, thickness-shear mode. Temperature level is 
maintained by a two-stage oven. The oscillator weighs approximately 



equipment manuals furnished by the manufacturers. I 1, 6, 
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three and one-half pounds; its dimensions are 4V x x 6 7/8" * 

The third type of oscillator tested, a Sulzer Laboratories 

Model 2.5 frequency standard, is also a solid-state unit. It has 

frequency outputs of two and one -half megacycles, one megacycle, and 

100 kilocycles. Stability as stated by the manufacturer is one part 
10 

in 10 or better per day. The crystal-unit is a two and one-half 
megacycle military-type designed by Warner of Bell Telephone Labora- 
< tories. It is mounted in a proportionally-controlled double oven. 

The oscillator weighs approximately six and one-half pounds; its 
dimensions are x 4V x 11 3/4". 
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3* Test Program. 

Purpose and Procedure 

The purpose of the testing program was to reveal the effects 
of static position, swinging motion, vibration, and shock on the 
frequency stability of the crystal oscillators described above. To 
accomplish this goal, each of the oscillators was subjected to certain 
conditions and combinations of position and motion while a constant 
observation was made of its frequency output. 

The frequency standard employed as a reference for the tests 
was a Varian Model V-4700 Rubidium Vapor Frequency Standard with a 
long term frequency stability of 5 x 10 ^in any one-year period 
(standard deviation). To insure the reliability of the standard 
during the test period, its output was continually phase -compared 
against the 20 kilocycle signal of WWVL. This comparison was made 
with a Textran Series 599 VLF Tracking Receiver. The stability 
during the period from 23 July 1965 to 1 November 1965 was 2*3 x 10 ^ 
(standard deviation). 

The output of the test oscillator was phase -compared against 
the Rubidium standard and the output of the phase comparator was 
recorded on a strip chart recorder. The basic test system is shown 
in Figure 1. Examination and processing of the recorded data 
produced the results presented in Section 5. 
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Circuits and Equipment 

While the test circuit was basically that shown in Figure 1, 
the laboratory hook-up was more complicated due to the simultaneous 
testing of as many as four oscillators. All phase comparisons were 
made at 100 kilocycles, necessitating the use of frequency dividers 
in the case of the Sulzer D5 ' s with their five-megacycle outputs. 

The phase comparators, frequency dividers, and oscillators also 
required separate power supplies. Figure 2 shows an extension of 
the basic test system for the case of four oscillators under test. 

A thermistor-bridge circuit and a chart recorder were employed 
to maintain a constant record of room temperature for the duration of 
the tests. 

Swinging motion was provided by a specially constructed swinging 
platform shown in Figures 3 and 4. The swing was driven by a three- 
quarter horsepower electric motor coupled to a hydraulic transmission. 
The transmission permitted adjustment of swing frequency and the 
slotted drive wheel permitted adjustment of swing amplitude. 

Vibration was accomplished with a Goodmans Model 390 A Vibration 
Generator. A low-frequency function generator and a power amplifier 
were used to supply the necessary input. The shake table with the 
three Sulzer D5 oscillators mounted on it is shown in Figure 5. 

Figure 6 illustrates the equipment used to provide a periodic 
lif t-and-release motion for the shock testing described in a later 
paragraph. The equipment consisted of a direct-current motor with 
a cam attached to its shaft, and geared for high-torque, slow-speed 
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FIGURE 1 

BASIC TEST SYSTEM 




FIGURE 2 

* 

BASIC TEST SYSTEM EXTENDED FOR 
SIMULTANEOUS TESTING OF FOUR SULZER OSCILLATORS 
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SWING PLATFORM, VIEW 1 




FIGURE 4 

SWING PLATFORM, VIEW 2 
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SULZER D5 OSCILLATORS MOUNTED ON SHAKE TABLE 
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FIGURE 6 

SHOCK TESTING EQUIPMENT 
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operation. The cam periodically lifted and released a hinged bar 
from which was suspended the mounting platform for the oscillators. 

Description of Tests 

Static Tests . The static tests consisted of operating the test 
oscillator for at least 48 hours on each of its six sides. The 
oscillator was operated first in its normal position, then rotated 
counterclockwise through the side, upside down, and side positions. 
It was then operated on each of the remaining two sides. After 48 
hours in each of these positions, the oscillator was returned to its 
normal position for an additional 48 hours of testing. These tests 
were conducted on all of the oscillators. Table I lists the various 
static tests and their descriptions. 

TABLE I 

STATIC TEST DESCRIPTIONS 



Test 



Oscillator Position 



Static One 
Static Two 
Static Three 
Static Four 
Static Five 
Static Six 



Oscillator in normal operating position 
Oscillator tipped onto its side 
Oscillator tipped upside down 
Oscillator tipped onto its side 
Oscillator tipped forward onto its face 
Oscillator tipped backward onto its back 
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Swing Tests, These tests were also conducted over a basic 



48 hour test period. The test oscillator or oscillators were 
fastened to the swing platform and the swing was then placed in 
motion. The Western Electric 0-76A/U was swung at various ampli- 
tudes and frequencies in one position only, the normal position, 
with the front panel of the oscillator parallel to the leading edge 
of the swing platform and perpendicular to the line of motion. The 
four Sulzer oscillators were swung simultaneously, and in two posi- 
tions, the first being identical to that described for the 0-76A/U. 
For the second position, the oscillators were rotated 90° in the 
plane of the platform. The period and amplitude of the various tests 
are given in Tables II and III. 



TABLE II 



SWING TEST DESCRIPTIONS, 
WESTERN ELECTRIC OSCILLATOR 



Test 



Amplitude* 



Period 



Swing Four 
Swing Five 
Swing Six 



Swing One 
Swing Two 



Swing Three 




30 seconds 
10 seconds 
30 seconds 
30 seconds 
10 seconds 
20 seconds 



*10°, 10° represents a symmetric swing of 10° off the vertical 
in each direction. 
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TABLE III 



SWING TEST DESCRIPTIONS, 
SULZER OSCILLATORS 



Test 


Amplitude* 


Period 


First Position: 






Swing One 


5° 10° 

n » o 


30 seconds 


Swing Two 


5 , 10° 


10 seconds 


Swing Three 


24°, 30° 


30 seconds 


Swing Four 


24 , 30 


10 seconds 


Second Position: 






Swing Five 


5°, 10° 


30 seconds 


Swing Six 


5°, 10° 


10 seconds 


Swing Seven 


24°, 30° 


30 seconds 


Swing Eight 


24 , 30 


10 seconds 



*5°, 10° represents an asymmetric swing of 5° off vertical in 
one direction and 10 off vertical in the opposite direction 



Vibration Test . This test was conducted only with the Sulzer 
D5 oscillators. The oscillators were mounted on a one-quarter inch 
thick aluminum plate with a diameter of twelve inches which was then 
attached to the vibration generator (see Figure 5). They were vi- 
brated simultaneously at a frequency of 20 cycles per second with a 
force of approximately one g for a period of 12 hours. After a 
12 hour quiet period, the oscillators were once again subjected to 
the same vibration condition for a period of ten hours. 

Shock Test . The shock, or impact, test was also conducted only 
with the Sulzer D5 oscillators. The test consisted of a free fall 
from a height of three-tenths of an inch onto asphalt tile laid on 
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concrete (see Figure 6). As in the vibration test, all three 
oscillators were tested simultaneously. Impacting occurred at a 
rate of once every ten seconds and was conducted continuously for 
12 hours. 

Environment . All tests were conducted in an open laboratory 
with room temperature normally ranging daily from 21°C to 23°C. Over 
the 43 day test period the mean temperature was 22°C with a low of 
18°C and a high of 29°C. 

Data furnished by Sulzer Laboratories with respect to effects 
of temperature change on frequency stability of the Sulzer oscillators 
is shown in Table IV. Similar data for the Western Electric oscil- 
lator was not available. 

TABLE IV 

EFFECTS OF TEMPERATURE ON 
FREQUENCY STABILITY OF SULZER OSCILLATORS* 



Oscillator 



Temperature 

Change 



Maximum Af* * 
F 



Model 2.5 

Frequency Standard 
Model D5 

Model D5 

Model D5 

Model D5 



ambient temperature + 2 x 10 

20°C changed + 15°C 

ambient temperature + 1 x 10 

-54°C vs + 25°C 



ambient temperature + 5 x 10 

-18°C vs + 25 °C 

ambient temperature + 5 x 10 

+60°C vs + 25°C 



ambient temperature + 1 x 10 

+65°C vs + 25°C 



^Information furnished by Sulzer Laboratories 

**Defined on Page 14 
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4. Data Acquisition and Processing. 

Acquisition 

During the course of each test, the outputs of the various phase 
comparators were fed to strip chart recorders which provided a 
permanent data record for processing and analysis at a later date. 
Temperature data was also permanently recorded. The phase -comparison 
data was in the form shown in Figure 7. 



/6©0 2oao 2400 o4a> 




SAMPLE OF RECORDED PHASE COMPARATOR OUTPUT 
The output of a phase comparator indicated the difference between 
the frequencies of the two input signals. If the frequencies of the 
signals were equal, the phase difference of the signals was constant, 
and the output of the phase comparator was constant. If the fre- 
quencies were not equal, however, the signals passed in and out of 
phase with each other periodically and the cyclic change' in phase 
difference produced a similar variation in the phase comparator output. 
By recording this cyclic output on a strip-chart moving at a regulated 
speed, the rate of phase shift, cycles/unit time, was determined. 

For example, as shown in Figure 7, during the period from 2100 to 2200 
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the phase comparator passed through 1.05 cycles; the rate of phase 
shift during this period was 1.05 cycles/hour. 

All phase comparison was done at a frequency of 100 kilocycles. 
Therefore, during the one -hour period, the standard signal went 

Q 

through 3600 x 100,000 or 3.6 x 10 cycles while the test signal 

8 

went through 3.6 x 10 - 1.05 cycles. 

The difference between the frequencies, expressed as a fraction, 

is called the fractional frequency error or the frequency offset. 

It is written as A f where 

F 

f *= frequency of test signal 
F = frequency of standard signal 
Af « f • F 

The frequency offset for the above example was 

A f = -1.05 = -2.9 x 10' 9 

F 3.6 x 10 S 

It is customary to refer to the amount of phase shift in terms 

of time rather than frequency. At 100 kilocycles, the time for one 

cycle is ten micro-seconds and 1.05 cycles of phase shift per hour 

is equivalent to 10.5 micro-seconds of phase shift per hour. The 

frequency offset may then also be computed as 

Af - A t = -10.5 x 10~ 6 = -2.9 x 10 -9 

F T 3600 

where T = period of observation in seconds 

£t = phase shift in seconds occurring during T 

Processing 

The recorded data for each test was examined in the manner out- 
lined above. The number of micro-seconds of phase shift per hour 
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was noted for each hour of each test period. These figures were 
used to compute hourly values of frequency offset, Af/F, which 
were then plotted versus time. A sample plot is shown in Figure 8. 

The positive values of frequency offset indicate that the frequency 
of the test oscillator was above that of the standard. The upward 
drift of the points represents an increase in frequency of the test 
oscillator. This is a normal condition, caused by aging of the 
crystal # 

Next, a digital computer was employed to fit a straight line 
to the plotted points and to measure the dispersion of the points 
about this fitted line by computing their standard deviation. The 
program used for these computations is given in the Appendix. The 
slope of the fitted line represented the frequency drift rate of the 
test oscillator and the standard deviation of the points was an 
indicator of the oscillator's frequency stability, with lower values 
of standard deviation indicating better stability. Figure 9 shows 
the sample plot presented in Figure 8 with a straight line fitted 
to the points. The choice of a one-hour interval between the points 
was purely arbitrary. Figures 10 and 11 present the same data plotted 
at intervals of one-half hour and two hours respectively. It may be 
noted that the smallest variation of Af/F plotted in Figure 9 
is 7 x 10 H ; in Figure 10, 14 x 10 ^ ; and in Figure 11, 3.5 x 10 
This is because the smallest value of Af/F is inversely propor- 
tional to the length of the observation interval. Using the 40-line 
chart paper shown in Figure 7, with the recorder calibrated for ten 
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micro-seconds of phase shift for full scale deflection, the smallest 
increment of phase shift readable is 0.25 micro-second. Therefore, 
for a one-hour observation interval, the smallest value of Af/F 
(other than zero) is 

Af = At = 0.25 x 10~ 6 = 7 x 10~ U 

F T 3600 

For the purposes of this investigation, the choice of a one-hour 
observation interval proved to be quite satisfactory. Additional 
accuracy could have been obtained, however, by using recording equip 
ment with 100-line paper, permitting phase shift to be read to 0.1 
micro-second. 

In summary, the data processing was accomplished in three steps 
First, the strip charts were examined for the amount of phase shift 
per hour. Second, this data was converted to values of frequency 
offset and plotted against time. Third, to the plotted points was 
fitted a straight line, whose slope was thus obtained, and a measure 
of the variation of the points about the line was computed. 
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5. 



Test Results. 



The results of all tests are presented in Tables V through 
IX and Figures 12 through 37. The complete testing program of 
static, swinging motion, vibration, and shock tests was conducted 
with only the three Sulzer Model D5 oscillators. The two frequency 
standards were not subjected to vibration and shock. 

Prior to the arrival of the Sulzer equipment, static and 
swinging motion tests were begun with two Western Electric fre- 
quency standards. One of these was the oscillator mentioned in 
Section 2, a model 0-76A/U, and the second was an earlier model, 
076/U. It was soon determined that the drift rate and frequency 
stability of the older 0-76/U were affected by fluctuations in room 
temperature. The combined weight of the two oscillators (approxi- 
mately 180 pounds) also proved to be too great a load for the swing 
platform and prevented simultaneous swing testing. A circuit 
failure in the 0-76/U, coupled with the problems of its poor stabi- 
lity and load on the swing, caused testing with this particular 
obcillator to be discontinued. The 0-76A/U, however, was unaffected 
by changes in room temperature and operated properly during the test 
period . 

The results of the static and swing tests for this oscillator 
are given in Table V and Figures 12 through 15. In some instances 
the movement of the oscillator from one static position to another 
created a temporary frequency disturbance and several hours were 
required for the drift rate to return to its former state. 
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During this recovery period the change in frequency offset was 
relatively rapid and therefore the total test period was divided 
into sections to permit better fitting of straight lines to the 
data points. The three highest values of drift rate plotted in 
Figure 12, and the three highest values of standard deviation plotted 
in Figure 13, occurred during the recovery periods. For reference 
purposes the oscillator was operated in the Static One or normal 
position at the beginning and end of both the static test series 
and swing test series. The M motor only" entry in Table V refers to 
a period when the swing drive was disconnected and the swing platform 
was at rest, but the drive motor was running. This particular test 
was conducted to check for possible effects of vibration generated by 
the motor and hydraulic transmission. 

After a three-week warm-up period, the new Sulzer equipment was 
tested. It was recognized that this was an inadequate warm-up time, 
but the late arrival of the equipment and the pressure of time made 
it mandatory to begin the tests. The results of these tests are given 
in Tables VI through IX and Figures 16 through 37. The room tempera- 
ture variations had no noticeable effect on the drift rate or fre- 
quency stability of these oscillators. 
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TABLE V 



TEST RESULTS FOR 

WESTERN ELECTRIC 0-76A/U OSCILLATOR 
SERIAL NO. 43 



Test 






Condition 


Length 

(Hours) 


Drift Rate „ 
Af/F (pp 10 ) 
hour 


Standard 
Deviation 
(x 10*) 


Static 1 






normal 


48 


0.110 


0.35 


Static 2 


(a) 




side 


24 


0.473 


1.05 


Static 2 


(b) 




side 


24 


0.141 


0.53 


Static 3 






upside down 


48 


0.143 


0.87 


Static 4 






side 


40 


0.172 


0.53 


Static 5 


(a) 




front 


12 


0.835 


1.12 


Static 5 


(b) 




front 


36 


0.053 


0.44 


Static 6 


(a) 




back 


6 


1.494 


1.11 


Static 6 


(b) 




back 


20 


0.237 


0.49 


Static 6 


(c) 




back 


22 


0.144 


0.43 


Static 1 






normal 


48 


0.102 


0.82 


Static 1 






normal 


192 


0.037 


0.57 


Static 1 






normal 


158 


0.025 


0.41 


Swing 1 


10°, 


10°, 30 


sec.* 


48 


0.133 


0.80 


Swing 2 


10 , 


10° 10 
n 7 


sec . 


48 


0.033 


0.80 


Swing 3 


20 °, 


20 , 30 


sec. 


48 


0.031 


0.72 


Swing 4 


5 °, 


5 o> 30 


sec . 


48 


-0.003 


0.86 


Swing 5 


5° 


5 a 10 


sec . 


48 


0.047 


1.29 


Swing 6 


10° , 


10 , 20 


sec. 


48 


-0.037 


0.81 


Static 1 


(motor only) 


normal 


48 


0.039 


0.67 


Static 1 




normal 


48 


0.127 


0.96 


Static 1 




normal 


48 


0.033 


0.60 


Static 1 




normal 


97 


0.025 


0.73 



*10 



o 

> 



the vertical 



10° , 30 sec. represents a symmetric swing of 10° off 
in each direction with a 30-second period. 



23 



tm 



:it: 



tit 

, n:r! 



ifH 

mi 



m 



ii'“ 






-t. 




lii - : 


m* 




-»-* ♦ 














, * 


* IT 1 


trt‘ 






tm 


! t tl 


ill' 


tip 


[ill 


: Itl! 


im 


lift f 
ffl-4 


ti tii titl tttt 

ii l|: f. if 


iHG 


mi 


n 






*xxi 


tfH 


ttti 


' t lit 


fH: 


* m 




|-pi; 


* i P : 


:i [ t 


HP 


im 


ntj 


!' ■['] 




tfu 


iin 

no 


.iiii 


i rn 


fjij 


i Hit 

_ i 1 1 


Ifl 




26 




TABLE VI 



TEST RESULTS FOR 

SULZER MODEL 2.5 FREQUENCY STANDARD 
SERIAL NO. 256 



Test 








Condition 


Length 

(Hours) 


Drift Rate .. 
Af/F (pp 10 l ) 
hour 


Standard 
Deviation 
(x 10“) 


Static 1 








normal 


65 


2.96 


5.04 


Static 2 








side 


49 


2.38 


6.18 


Static 3 








upside down 


49 


2.05 


5.41 


Static 4 








side 


48 


1.79 


4.17 


Static 5 








front 


48 


1.68 


5.29 


Static 6 








back 


48 


1.53 


11.35 


Static 1 








normal 


68 


1.47 


9.30 


Swing 1 


5 o> 


10°, 


30 


sec .* 


49 


1.26 


7.95 


Swing 2 


5 , 


10°, 


10 


sec. 


50 


1.16 


9.40 


Swing 3 


24 0> 


30° , 


30 


sec. 


47 


1.04 


9.07 


Swing 4 


24, 


30 , 


10 


sec. 


47 


1.09 


6.97 


Swing 5 


5 o’ 


10°, 


30 


sec . 


48 


0.91 


3.91 


Swing 6 


5 o> 


10°, 


10 


sec. 


48 


0.98 


3.81 


Swing 7 


24, 


30 , 


30 


sec. 


48 


0.70 


4.46 


Swing 8 


24°, 


30 , 


10 


sec . 


48 


0.77 


4.83 


Static 1 








normal 


48 


0.81 


3.76 



*5°, 10°, 30 sec. represents an asymmetric swing of 5° off 
vertical in one direction and 10 off vertical in the opposite 
direction with a period of 30 seconds. 
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TABLE VII 



TEST RESULTS FOR 
SULZER MODEL D5 OSCILLATOR 
SERIAL NO. 25 



Test 




Condition 


Length 

(hours) 


Drift Rate . . 
Af/F (PP 10 ) 

hour 


Standard 
Deviation 
(X 10 ir ) 


Static 1 




normal 


65 


-4.14 


8.46 


Static 2 




side 


49 


-3.61 


7.28 


Static 3 




upside down 


49 


-3.13 


6.17 


Static 4 




side 


48 


-2.81 


7.10 


Static 5 




front 


48 


-2.51 


5.93 


Static 6 




back 


48 


-2.26 


8.22 


Static 1 




normal 


68 


-2.31 


5.66 


Swing 1 


5<\ 


10°, 30 sec.* 


49 


-2.73 


15.06 


Swing 2 


5°, 


10°, 10 sec. 


50 


-2.25 


12.63 


Swing 3 


24° 


30 , 30 sec. 


47 


-1.90 


12.63 


Swing 4 


24 , 


30°, 10 sec. 


47 


-1.47 


4.28 


Swing 5 


5°, 


10°, 30 sec. 


48 


-1.41 


6.64 


Swing 6 


5°, 


10° , 10 sec. 


48 


-1.42 


9.28 


Swing 7 


24 , 


30 , 30 sec. 


48 


-1.61 


11.23 


Swing 8 


24°, 


30°, 10 sec. 


48 


-0.68 


9.30 


Static 1 




normal 


48 


-0.81 


10.22 


Vibrate 1 




20 cps, 1 g 


12 


-5 .02 


3.92 


Static 1 




normal 


12 


-3.59 


7.26 


Vibrate 2 




20 cps, 1 g 


10 


-2.42 


9.64 


Static 1 




normal 


18 


-3.01 


8.45 


Static 1 




normal 


17 


-1.01 


7.99 


Shock 1 




0.3" fall 


12 


+0.71 


9.11 


Static 1 




normal 


60 


-1.63 


4.89 



*5°, 10°, 30 sec. represents an asymmetric swing of 5° off 
vertical in one direction and 10 off vertical in the opposite 
direction with a period of 30 seconds. 
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TABLE VIII 



TEST RESULTS FOR 
SULZER MODEL D5 OSCILLATOR 
SERIAL NO. 19 



Test 


Condition 


Length 

(hours) 


Drift Rate.. 
Af/F (pp 10 ) 

hour 


Standard 
Deviation 
(x 10 L1 ) 


Static 1 


normal 


65 


-6.29 


8 , 46 


Static 2 


side 


49 


-4.82 


8.33 


Static 3 


upside down 


49 


-3.96 


8.36 


Static 4 


side 


48 


-3.52 


8.32 


Static 5 


front 


48 


-3.02 


8.86 


Static 6 


back 


48 


-3.28 


8.98 


Static 1 


normal 


68 


-2.86 


6.32 


Swing 1 


5°, 10°, 30 sec.* 

A 7 A 7 


49 


-3.31 


8.75 


Swing 2 


5 , 10 , 10 sec. 

o * o 5 


50 


-2.68 


7.68 


Swing 3 


24 , 30 , 30 sec. 

n 7 n J 


47 


-2.46 


4.74 


Swing 4 


24 , 30 , 10 sec. 

n 7 n 7 


47 


-2.22 


4.12 


Swing 5 


5 , 10 , 30 sec. 

o 7 o 7 


48 


-2.03 


6.06 


Swing 6 


5 , 10 , 10 sec. 

n 7 n 7 


48 


-2.12 


4.67 


Swing 7 


24 , 30 , 30 sec. 

O 7 A 7 


48 


-2.12 


5.88 


Swing 8 


24 , 30 , 10 sec. 


48 


-1.77 


9.96 


Static 1 


normal 


48 


-1.64 


8.22 


Vibrate 1 


20 cps, 1 g 


12 


-4.63 


9.55 


Static 1 


normal 


12 


-3.11 


4.25 


Vibrate 2 


20 cps, 1 g 


10 


-3.81 


5.10 


Static 1 


normal 


18 


-3.03 


9.52 


Static 1 


normal 


17 


-1.69 


5.86 


Shock 1 


0 . 3" fall 


12 


-2.36 


8 .45 


Static 1 


normal 


60 


-2.19 


10.06 


*5° 


, 10°, 30 sec. represents an 


asymmetric swing 


of 5° off 


vertical 


in one direction and 


10 off 


vertical in the 


opposite 


direction 


with a period of 30 


seconds . 
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TABLE IX 



TEST RESULTS FOR 
SULZER MODEL D5 OSCILLATOR 
SERIAL NO. 22 



Test Condition Length Drift Ra£| Standard 

(hours) Af/F (pp 10 ) Deviation 

hour (x 10 ) 

— 



Static 1 




normal 


65 


-6.68 


4.97 


Static 2 




side 


49 


-5.14 


7.44 


Static 3 




upside down 


49 


-4.31 


6.31 


Static 4 




side 


48 


-4.26 


10.12 


Static 5 




front 


48 


-2.90 


9.77 


Static 6 




back 


48 


-2.64 


7.67 


Static 1 




normal 


68 


-2.26 


8.64 


Swing 1 
Swing 2 


5°, 

5° 

A ? 


10°, 30 sec. * 
10°, 10 sec. 


49 

50 


-2.41 

-2.10 


7.85 

5.20 


Swing 3 


24 , 


30 , 30 sec. 


47 


-1.76 


4.09 


Swing 4 


24°, 


30°, 10 sec. 


47 


-1.54 


4.08 


Swing 5 


5°, 


10 , 30 sec. 


48 


-1.42 


4.47 


Swing 6 


5°, 


10°, 10 sec. 


48 


-1.32 


3.07 


Swing 7 


24°, 


30°, 30 sec. 


48 


-1.39 


3.09 


Swing 8 


24° , 


30°, 10 sec. 


48 


-1.08 


4.65 


Static 1 




normal 


48 


-0.95 


5.00 


Vibrate 1 




20 cps, 1 g 


12 


-4.50 


7.46 


Static 1 




normal 


12 


-2.74 


2.18 


Vibrate 2 




20 cps, 1 g 


10 


-2.80 


2.56 


Static 1 




normal 


18 


-1.99 


4.26 


Static 1 




normal 


17 


-0.55 


3.32 


Shock 1 




0.3" fall 


12 


-0.02 


7.39 


Static 1 




normal 


60 


-1.14 


4.37 



*5°, 10°, 30 sec. represents an asymmetric swing of 5 off 
vertical in one direction and 10 off vertical in the opposite 
direction with a period of 30 seconds. 
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1 
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6. Observations. 



At the onset of the investigation, the investigator had hoped 
to be able to correlate any observed changes in drift rate or frequency 
stability with the conditions of the particular test during which the 
change was noted. A seemingly ideal situation was present with the 
four Sulzer oscillators, three being identical models and the fourth 
an excellent frequency standard. It had been planned to obtain addi- 
tional information from Sulzer Laboratories concerning the orienta- 
tion of the crystal plane within the oscillators. With this informa- 
tion the investigator had hoped to relate the results of the static 
tests with crystal position. 

The manufacturer stated, however, that the position of the crystal 
plane is not fixed, but varies with each individual oscillator. 

Since the only remaining method of determining the crystal orienta- 
tion would have required disturbing the double-oven assemblies, this 
second aspect of the investigation was abandoned. 

The correlation of changes in frequency stability and drift rates 
with test conditions in the case of the Sulzer equipment was hampered 
by the fact that the limited time available for testing required that 
the tests be conducted during the initial warm-up and aging period 
of the oscillators. As the equipment was new, no data record was 
available for use as a behavior reference. These conditions prevented 
the investigator from making the detailed conclusions originally hoped 
for. Instead, general observations of the behavior of the oscillators 
while under test were made. These observations are: 
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Western Electric Frequency Standard 0-76A/U 

9 

Static Tests . Drift rates varied between 0.04 (pp 10 ) per 
9 

hour and 1.49 (pp 10 ) per hour, with the lowest value occurring 
while the oscillator was in the normal position. However, each time 
the oscillator was in this position, a different drift rate was 
observed (0.11, 0.10, 0.04). During a portion of the Static Five 
test with the oscillator tipped forward on its face the drift rate 
was 0.05 (pp 10 9 )/hr. 

-9 

Stability (standard deviation) varied from 0.35 x 10 to 
-9 

1.12 x 10 with the lowest value occurring while the unit was in a 
normal position. However, as in the case of the drift rates, each 
time the oscillator was returned to its normal position, a different 
value was observed (0.35, 0.32, 0.57). Once again there were instances 
where the stability while the oscillator was in some position other 
than normal was better than that observed for some of the normal 
position test periods. 

The maximum values of drift rate and standard deviation occurred 

during the recovery periods as discussed in Section 5. 

Swing Tests . The drift rate was certainly affected by the 

swinging motion although during the test period from 12 to 14 

September while the oscillator was in the Static One position the 

drift rate was almost as great as the largest observed during 

9 

swinging. The rates varied from 0.133 (pp 10 )/hr. to -0.037 
(pp 10 9 )/hr. 
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Values of standard deviation were generally greater during 
swinging than while at rest, indicating poorer stability under these 
conditions . 

Sulzer Model 2 . 5 Frequency Standard 

Static Tests . The drift rate gradually decreased as the crystal 
aged and no apparent effects of static position were observed c 

Stability varied between 4.2 x 10 ^ and 11 A x 10 ^ with the 
greatest value observed during the Static Six test, while the unit 
was tipped on its back. 

Swing Tests . The drift rate continued to decrease gradually 
with further aging of the crystal. In three out of four instances, 
however, the drift rate increased slightly when the swing frequency 
was increased. 

Stability improved as the swing tests progressed, with values 
of the standard deviation decreasing from 9.3 x 10 ^ to 3.8 x 10 
In three out of four instances, however, the standard deviation in- 
creased when the amplitude of the swing was increased. 

Sulzer Model D5 Oscillators 

Static Tests . As in the case of the Sulzer frequency standard, 
the drift rate for each of the three Model D5 oscillators gradually 
decreased during the aging period. 

-11 -11 

Stabilities varied between 5 x 10 and 10 x 10 with no 
common behavior pattern evident. One oscillator exhibited a downward 
trend, one an upward trend, and one remained fairly constant. 
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Swing Tests , The drift rates generally decreased during the 



swinging tests as aging continued. 

-11 -11 

Stabilities varied from 3 x 10 to 15 x 10 . The lowest 

values of standard deviation for each oscillator were obtained during 
one or another of the swing tests and not during a static period. 

In three out of four instances when the swing frequency was increased, 
the stability of various oscillators improved. The exceptional cases, 
however, did not occur simultaneously. 

The stability of two of the oscillators, serial no. 19 and 
serial no. 22, improved in three out of four instances when the 
amplitude of the swing was increased. 

Vibration and Shock Tests . All drift rates were increased by 
three or four parts per 10^ per hour during the initial vibration 
test. 

Standard deviation varied between 2.2 x 10 ^ and 10.1 x 10 
All standard deviations were increased during the shock test; the 
amounts of increase varied from one to four parts per 10^\ 

In conclusion, the drift rate and frequency stability of the 
frequency standards were apparently affected to a small degree by 
swinging motion, although this motion had no definite effect on the 
Sulzer Model D5 oscillators. The D5 oscillators were affected 
slightly by shock and vibration. In only a few cases, however, such 
as the change in stability during the shock test, or the change in 
drift rate during the initial vibration test, did all three Model D5 
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oscillators respond in a similar manner. Generally, opposite re- 
actions occurred, i.e., during a given test period, drift rates of 
two oscillators would increase while the third decreased, or the 
stabilities of a different pair of oscillators would decrease while 
the third increased or remained constant. 

In spite of the lack of positive correlation of changes in 
drift rate and stability with conditions of position and motion, 
the results of the investigation were fruitful in that they indicated 
a lack of vulnerability of the oscillator units to motion in general. 
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7. Recommendations for further study. 

Since the initial warm-up and aging period of the oscillators 
is now over, and a backlog of behavior data has been accumulated, 
further investigation should produce more positive results. A repe- 
tition of the tests described in Section 3 would be in order, as 
would be further testing with the Sulzer D5 oscillators mounted in 
shock mounts. 

Test periods should be at least 48 hours in length and, if 
time permits, each test should be preceded and followed by a similar 
amount of time with the oscillators in a normal position. Record- 
ing equipment utilizing 100-division paper should be used to record 
the phase comparator output so that with ten micro-seconds full- 
scale deflection, phase shift at 100 kilocycles may be read to one- 
tenth or one-twentieth of a micro-second. 
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APPENDIX 



A COMPUTER PROGRAM FOR FITTING A STRAIGHT LINE 
TO DATA POINTS AND COMPUTING THE STANDARD DEVIATION 
OF THE POINTS ABOUT THE FITTED LINE 
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oooooooooooooooooooooo 



PROGRAM LINEFIT 

THIS PROGRAM FINDS THE SLOPE AND Y-INTERCEPT OF THE BEST STRAIGHT 
LINE (IN THE LEAST SQUARES SENSE) FITTED TO A SET OF DATA POINTS. 
THE X DATA ARE ASSUMED TO BE WITHOUT ERROR. TRANSFORMATIONS TO 
EITHER OR BOTH VARIABLES CAN BE MADE. THE PROGRAM WILL ALSO 
COMPUTE THE STANDARD DEVIATION OF THE POINTS. 

USAGE 

1. calling sequence 

MAIN PROGRAM. THE FOLLOWING DATA CARDS ARE READ. 

A. FORMAT (314) - FIRST CARD 
N * RIGHT-JUSTIFIED NUMBER OF: DATA POINTS. 

ICODEX,. TRANSFORMATION CODES FOR X 
ICODEY, TRANSFORMATION CODES FOR Y 

B. FORMAT (12F6.0) - SUBSEQUENT ' CARDS . 

X(1)»Y(1)*X(2)»Y(2)» ETC. 

C. TRANSFORMATION CODES 
=1 NONE 

=2 VARIABLE REPLACED BY SQUARE) ROOT OF VARIABLE 
=3 VARIABLE REPLACED BY SQUARE: OF VARIABLE 
= 4 VARIABLE REPLACED BY LOG TO BASE 10 OF VARIABLE 
=5 VARIABLE REPLACED BY LOG TO BASE E OF VARIABLE 

11111 

TEST PERIOD FROM 1700 4 OCT TO 1300 7 OCT, SI, STATIC ONE, FINAL 
DIMENSION X( 200 ) ,XSAVE(200 ) , Y(200 > , YSAVE(200 ) , YF I T ( 20 0 ) , D I FFS ( 20 0 
104 READ 999* N, ICODEX, ICODEY 
999 FORMAT ( 3 I 4 ) 

PRINT 998, N, ICODEX, ICODEY 

998 FORMAT ( 4H1N = I 4 , 3X8H I CODEX = I 2 , 3X8H I CODE Y =12//) \ 

READ 100, (X( I),Y(I ), ‘ \ 

100 FORMAT (12F6.0) 



XSA VE ( I ) =X ( I ) \ \ 

996 YSAVE ( I ) = Y ( I ) 

YSUM =0. * \ 

SUMX = 0 • 

SUMY = 0. 

SUMX Y = 0. 

SUMXX = 0. 

CALL TRANS(N,X,Y, ICODEX, ICODEY) 

XN = N 

DO 30 I M, N 
SUMX=SUMX+X( I ) 

SUMY = SUMY + Y ( I ) 

SUMXX = SUMXX ♦ ( X ( I ) *X ( I ) ) 

30 SUMXY=SUMXY+ <X(1)*Y(I>) 

DENOM = XN*SUMXX-SUMX*SUMX 

ANUM = SUMY * SUMXX - SUMX * SUMXY 

A = ANUM / DENOM 

BNUM= XN*SUMXY-SUMX*SUMY 

B = BNUM / DENOM 

DO 40 J = 1 , N 

YC = A + B * X(J) 

• DIFF = Y( J) - YC 
4o YSUM = YSUM + D I FF*D I FF 

S2 = YSUM / (XN - 1 . ) 

S I G=S2/DENOM 
SRA=SQRTF(SUMXX*SIG) 

SRB=SQRTF(XN*SIG) 

PRINT 997, A, B, SRA , SRB 

FORMAT ( 24H0 STRAIGHT L I NE’ F I TT I NG// 25H USING THE FORM Y*BX «• A, 
1/ 4H A = E14.8/ 4H B = E14.8/ 26H STANDARD DEVIATION OF A = E14.8 



DO 996 I=i, N 
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997 



2/ 26H STANDARD DEVIATION OF B = E14.0) 

SUMF I T = . 0 
SD I FFS=0 . 

S D 1 FFSQ = 0 . 

DO 90 1=1, N 

YF I T ( I ) =B*X ( I ) ♦ A 

DIFFS( I )=YFIT( I )-Y( I ) 

SDIFFS = SDIFFS + DIFFS( 1 ) 

SDIFFSO^SDIFFSO + DIFFSCI )**2 
9o SUMF I T=SUMF I T+ DIFFS(I) 

995 PRINT 99»SUMFIT, (X( I ),YSAVE( I )>Y( I ),YFIT ( I ) # D I FFS < I ), 1 = 1, N) 

99 FORMAT (9H0SUMF IT =E12 ,5///10XlHX14XlHY9X10HYTRANSPOSE8X4HYF I T11X5. 
1DIFFS///(5E15.5) ). 

SIGMA = SQRTF((SDIFFSQ-SDIFFS**2/XN)/(XN -1:.)) : 

PRINT 1000, SIGMA 

1000 FORMAT (1H0,7HSIGMA =E15.8) 

GO TO lo4 

106 STOP l ■ 

END 

SUBROUTINE TRANS ( N , X, Y, ICODEX, ICODEY) 

DIMENSION X ( 200 ) , Y(200> 

GO TO ( lo# 20, 3o, 4 o # 5o>» ICODEX 
20 DO 21 I = 1 , N 
21. X( I ) =SQRTF ( ABSF ( X ( I ) ) ) 

GO TO 10 

30 DO 31 I = 1,N ' 

31 X(I) = X(I> * X(I) 

GO' TO 10 

40 DO 4l I = 1 , N 

41 X( I )* = LOG10F(X(1) ) • 

GO TO 10 

50 DO 5l I = 1 , N 

51 XII ) =LOGF ( X ( I ) ) 

10 GO TO (100, 200, 300, 400, 500), ICODEY 

200 DO 201 I = 1 , N 

201 Y(I)=SQRTF(ABSF(Y(I))> 

100 RETURN 

300 DO 3ol I = 1 , N 

301 Y ( I ) = Y(I) * Yd) 

RETURN 

400 DO 4 0 1 I = 1 , N 

401 Yd )=LOG10F( Yd ) ) 

RETURN 

500 DO 5ol I = 1 , N 



501 Y( I )=L06F(Y( I > ) 
END 
END 
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N 



68 ICODEX = 1 ICODEY = 1 



straight line, fitting ' 

USING THE FORM Y = 8X + A, 

A =-.9o372827E+o3 
B =- . 23o75352E+01 

STANDARD DEVIATION OF A = . 13877674E + 01 
STANDARD DEVIATION OF 8 = .349629Q2E-01 



SUMFIT = .10431E-06 



X 



•lOOOOE+oi 
•20000E^01 
• 3 0000E + 01 
• ^0 0 0 0 E + 0 1 
• 5 0000E^01 
.60000E+01 
. 70000E+01 
•80000E+01 
. 90000E+01 
•l0000E + 02 
.11000E+02 
.12000E+02 
. 13000E+02 
. 14000E+02 
. 15000E+02 
•16000E+02 
. l7000E"-02 
. 18000E+02 
. 19000E+02 
•20000E+02 
• 21000E+02 
.22000E+02 
.23000E+02 
. 24000E+02 
.25000E+02 
. 26000E+02 
. 27000E+02 
. 28000E+02 
. 29000E+02 . 
.30000E+02 
•31000E+02 
.32000E+02 
.33000E+02 
. 34000E + 02 
. 35000E+02 
.36000E+02 
•37000E+02 
. 38000E + 02 
. 39000E‘*-02 

• A0000E + 02 
.41000E+02 

• A2000E^02 



Y 



- . ^o^ooE+O 3 

“ . 9i7qoE+o3 

- . 91000E + 03 

- .9l000E-*-03 

- .924qoE + o3 

- .9l700E-*-03 

- .91000E + 03 
- .93100E+03 
- . 93100E+03 
- . 93100E+03 

- .93100E + 03 
- . 93100E+03 
- . 93100E+03 
~ . 93lo OE-t-03 
- . 938ooE+o3 
- . 944o 0E + 03 
- . 94400E+03 
“ . 944q OE + O 3 
- . 94400E+03 
- . 944qoE+o3 
- . 95800E+03 
- . 965o 0E + 03 
- . 965o 0E + 03 
- ,958ooE+o3 
- . 965ooE+o3 
- . 9580QE+03 
- . 96500E+03 
- . 96500E+03 
~.972oOE+03 
.-.965ooE + o3 

- .97900E + 03 
-.97200E+03 
-.97900E+03 

- .965obE'4-03 

- .97900E + 03 

- .99300E + 03 
-.98600E+03 
- .99300E+03 
-.98600E+03 
- . 993o 0E + 03 
-.993oOE+o3 

- .993qoE + o3 



ytranspose 



- . 9 0 4 0 0 E «• 0 3 
-,9l7 0 oE+o3 
" .91000E+0 3 
- . 91 0 0 0 E^O 3 
-.924 0 oE+o3 
- . 91700E+03 
- .91000E*0 3 
-.93looE+o3 
~ . 93100E+0 3 
" . 9 3 1 o 0 E + 0 3 - 
- . 931o 0 E + 0 3 
- . 93100E + 03' 
" . 93lo OE + O 3 
- . 93lo 0E + 0 3 
- . 938o 0 E + 0 3 
" . 944ooE+o3 
- . 944qoE+03 
- . 944qoE+o3 
-.944 00 E + 0 3. 
-,944 00 E+o3 
- . 958ooE + o 3 
“ . 965 ooE + o3 

- . 965ooE + o3 
• . 958 Q oE+o3 
-.965qoE+o3 

95800E+0 3 
- , 965ooE+o3 
- . 965ooE+o3 
- . 972o 0 E + 0 3 
- . 96500E+03 
-.97900E+03* 
••.97200E + 03 
-.97900E+03 
- . 965o 0E + 03 
"#97900E+03 
-.99300E+03- 
-.98600E+03 
- . 993o 0 E + 03 
-.986ooE+o3 
-t993ooE+o3 
- . 993o 0 E + 0 3 
" .993ooE*o3 
JA 



Y F I T 



- . 9o6o4E + o3 
- . 9o634E + o 3 
• -.9l 0 65E+ 0 3. 

-.91296E+03 
-.91527E+03 
- . 91757E + 03 
- . 91988E + 03 
- . 92219E + 0 3 
-,9245oE+03 
- . 92680E + 03 
-.92911E+03 
- # 93142E + 03 
-.93373E+03 
- . 936Q3E+03 
- . 93834E+03 
• . 94o65E+o3 
-.94296E+03 
-.94526E+03 
-.94757E+03 
-.94988E+03 
- . 95219E+03 
-.95449E+03 
- . 95680E+03 
-.95911E+03 
-..96142E + 03 
-.96372E+03 
- . 966 o3E + o3 
- . 96834E + 03 
-.97065E+03 
-.97295E+03 
- . 97526E+03 
- , 97757E+03 
- . 97988E + 0 3 
- . 98218E+03 
- . 98449E+03 
-.98680E+03 
~ . 98911E+03 
• . 99141E+03 
-.99372E+03 
- . 996Q3E+03 
- . 99834E + 0 3 
-•100Q6E+04 



DIFFS 



- , 2o358E + ol 
.86567E+01 
~.65o88E+oO 
-.29584E+01 
,8734 ie + 0 1 
-.57348E+00 
~<988l0E+0l 
.88114E+01 
.65039E+01 
.41964E+01 

• 18888E + 01 
-.41869E+00 
-.27262E+01 
- . 5o338E+Ol 
-,34l3oE+00 

..33512E + 01 
. 10436E + 01 
- . 12639E + 01 
" . 357l4£+ol 
-.58790E+01 
.58135E+01 

♦l0 5 0 6 E+02 
.81984E+01 
-.11091E+01 
, 35833E + 01 
-.57242E+01 
-.10317E+01 
-.33393E+01 
.13532E+01 
-,79543E+0l 

• 37381E + 01 
-.55694E+01 
- . 87693E+00 
-.17184E+02 
-.54920E+01 

. 62005E+01 
- .31071E + 01 
♦15854E+01 
- • 77221E+01 
-.30297E+01 
-.53372E+01 
-.76447E+01 



. 43000E+02 
. 44q00E + 02 
.45000E+02 
. 46000E+02 
. 47oOOE^o2 
. 48000E + 02 
. 49000E*02 
. 50000E+02 
.51000E+02 
. 52000E*02 
• 53000E+02 
.54000E*02 
. 55000E+02 
.56000E+02 
.57000E+02 
. 58000E*02 

• 59oOOE-*-o2 
• 60000E+02 
• 61000 E +02 

• 6200 0E + 02 
. 63000E+02 
.64qooE*o 2 
.65oOOE + o2 
. 66 q00E + 02 
.67q00E+02 
. 60QOOE+O2 



- . 10000E+04 
- .lOOOOE+O 4 

- .10070E + 0 4 
- .l01 4 0E+0 4 

- .10070E + 0 4 

- .l0l 4 0E + 0 4 
- .10280E+0 4 
-.10280E+0 4 
- .10210E+04 

- .10280E + 04 
- . 10210E+0 4 
-.10210E+04 
-.10350E+0 4 
-.10350E+04 
- .10420E+04 

- .10420E + 04 
r . lo 4 2oE«-0 4 
- . l0 4 20E+0 4 
- . l0 4 20E+0 4 
* .10 4 2oE«-0 4 
- .10 4 2oE+0 4 

- .l0 4 90E + 0 4 
- .l0 56 0E+0 4 
- . 10 6 30E+0 4 
- .10690E+0 4 
- .10 6 3 qE*0 4 



- . 10000E+04 
- .10000E+0 4 ’ 

- . 1007oE^0 4 

- . 10l 4 0E+0 4 

- .10070E + 0 4 
- ,10l 4 0E+0 4 
-.10280E-0 4 

- . l0280E-*-0 4 
- . 10210E+0 4 
- .1O20OE+O 4 
-.10210E*0 4 

- .10210E + 0 4 ' 

- , 10350E + 0 4 

- .10350E*0 4 

- ♦ 1 0 4 2o E + 0 4 ' 

- . 10 4 2oE + 0 4 
- . 1o 4 2oE + o 4 
- . 10 4 2oE + 0 4 
- . 10 4 2oE«-o 4 
- . 10 4 2oE-*-o 4 

- .10 4 2oE + 0 4 

- ,10 4 9oE*0 4 

-.l0 56 0E+0 4 

- .10 63 oE«-0 4 ' 
-,10 6 9oE+0 4 

”«10^3qE+q4 



“tl0030E+04 

- .100 5 3E + o 4 

- . 10076E + 0 4 
-.10099Et-0 4 

- .10122E + 0 4 

- . 101 4 5E + 0 4 
-.10l68Ef04 

- .10191E + 0 4 

- .10214E-t-04 

• 10237E + 04 

7 . 10260E + 0 4 
f . 1O203E+O 4 
- . 10306E+0 4 

- . 10330E + 04 

- . 10353E + 04 
”.10376E-*-0 4 
- . 10399E + 0 4 

- .10 4 22E + 0 4 

- . lo 44 5E-*-o 4 
~ . 1 0 4 68E + 0 4 

- .10 4 9lE + 0 4 

- .10514E + 0 4 

- .l0 537 E + 0 4 
- . 1056 qE + 0 4 ' 
-.10583E+0 4 
-. 106 q 6E+0 4 



-.29523E+01 
-.52598E+01 
- , 56736E+00 
.41251E+01 
-,5l824E*0l 
- , 40996E+OO 
,112o3e*o2 
.88950E+01 
-, 4 l257E+00 
,42799E*0l 
-,5o276E*Ol 
- , 73352E+01 
, 43573E+0 1 
,20 4 98E+0l 
,67'j}22E + 0l 
.44?47E+01 
,2l?72E*0l 
-.1&038E+00 
- . £4®79E + 0 1 
-♦47955E+01 
7.1030E + 01 
-■.2|1o5e* 01 
• 2<J8i9E + 0l 
, 69744E+01 
, 1 o667£ + 0 2 
.23593E+01 



SIGMA = .5650916OE+O1 
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